Preprint typeset in JHEP style - HYPER VERSION 



SU-ITP-12/35 



Multiverse rate equation including bubble collisions 



Michael P. Salem 

Stanford Institute for Theoretical Physics and Department of Physics, Stanford University, 
Stanford, CA 94305, USA 



Abstract: The volume fractions of vacua in an eternally inflating multiverse are described by 
a coarse-grain rate equation, which accounts for volume expansion and vacuum transitions via 
bubble formation. We generalize the rate equation to account for bubble collisions, including 
the possibility of classical transitions. Classical transitions increase the connectedness of the 
landscape — with possible implications for the so-called staggering issue — however the basic 
intuition developed in prior literature for understanding the aforementioned volume fractions 
is easily updated to account for the effects. Whether or not our vacuum is likely to have been 
established by a classical transition depends on the detailed relationships among transition 
rates in the landscape. 
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1. Introduction 

We might live in an eternally inflating multiverse. Eternal inflation occurs whenever a suffi- 
ciently large volume is in a state sufficiently close to a vacuum in which the energy density is 
positive and the decay rate is smaller than the Hubble rate 0] , and /or whenever a sufficiently 
homogeneous field configuration evolves in a sufficiently flat, positive interaction potential 
[||, [|. These statements contain a number of qualifications, so it is worth noting that the 
first set of conditions is satisfied by the observed state of the local universe, if the so-called 
dark energy is due to vacuum energy, while the second set is similar to the initial conditions 
implicit in the simplest models of slow-roll inflation, though with the inflaton further up its 
potential Q. When it occurs, eternal inflation generates an endless spacetime in which every 
phase of vacuum takes place in a fractal mosaic of widely separated domains [||, || , the various 
vacua being attained either by bubble formation [^, ||, by stochastic diffusion Q, and/or by 



other processes [K], 11, [Eg ]. 

In an eternally inflating multiverse, the fraction of the global spacetime volume occupied 
by any of the various vacua is not directly observable. Nevertheless, these volume fractions 
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might be relevant to understanding the local conditions in our universe. For instance, the 
proper theory of initial conditions might be a theory of the multiverse as a whole, with our 
local "initial" conditions — i.e. the conditions describing the onset of slow-roll inflation in the 
particular phase of vacuum that gives rise to our universe — being selected according to their 
prevalence in the global spacetime. It is also possible that these volume fractions express a 
holographic dual to bulk spacetime physics |[3], [14], [15], [l(| 17, |T§| ]. 



In a multiverse where vacuum transitions occur predominantly via bubble formation, the 



volume fractions occupied by the various vacua are described by a rate equation || 19 [. An 
important caveat is that the predictions of the rate equation depend on the choice of global 
time foliation. This is because different global time foliations explore different regions in the 
diverging spacetime at different rates, while the exponential expansion of eternal inflation 
ensures that most of the total volume is near the boundary at any finite time cutoff. This 
expresses the measure problem of eternal inflation (for some recent reviews, see for example 



[ PQ , 21, |23|) Resolving the measure problem is of fundamental importance to (eternal) 
inflationary cosmology, but it is tangential to the thrust of this paper. 

This paper concerns another shortcoming of the rate equation, which is that it ignores 
bubble collisions. Semi-classical vacuum transition rates are exponentially suppressed, and a 
contribution to the rate equation from a bubble collision should involve a product of two such 
rates (one for each bubble in the collision), and so one might argue that bubble collisions can 
be ignored by expanding in powers of transition rates. However, these rates are exponentially 
staggered, meaning the product of two transition rates could be much larger than another 
transition rate. Moreover, over the course of its evolution each bubble collides with a diverging 
number of other bubbles. This divergence is regulated by the aforementioned measure, but 
it is not a priori clear how this resolution will play out. 

We perform our analysis with reference to both the scale- factor cutoff measure [24, 25, 
27] and the lightcone-time cutoff measure [^, 16]. Although the detailed effects of bub- 
ble collisions on vacuum volume fractions might be relevant to the aforementioned efforts to 
construct holographic duals to bulk dynamics, our analysis focuses on the qualitative effects 
of collisions on multiverse phenomenology. As such, we are uninterested in order-unity mod- 
ifications to the volume fractions. Instead, we focus on the effects of "classical transitions," 
whereby the kinematics of a bubble collision cause a transition to a new vacuum that might 



otherwise be relatively difficult to reach [11, |29| , 



For those very familiar with the literature (including the standard notation and assump- 
tions), our conclusions are simple to state. After including the effects of bubble collisions, the 
rate-equation transition matrix becomes 



Mij = Kij - 6ijKi + Jiklj^ijKkj , (1.1) 



kJ 



where the "fikij are related to the average volume fraction of vacuum i in the causal futures 
of collisions between bubbles of vacua k and £ which nucleate in vacuum j. The effects of 
bubble collisions are most significant when the first two terms in (|1.1| ) are zero but (because 
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of classical transitions) the third term is not. Meanwhile, the components of the dominant 
eigenvector Sj of My are to, leading order, 

E K «Pi + Ylj,k^ijkpi K jpi K kp! Kp n l + zJj.fe lp n jklKjlKkl 
4 X...X * , 1.2) 
Ki-q K Pn -q 

{Pa} 

where the sum covers the sequences of transitions that connect the dominant vacuum "1" to 
the vacuum i using the fewest number of upward transitions. Again, the effects of classical 
transitions can be significant because the first term in each factor can be zero when the second 
term is not. Note that the dominant vacuum — defined as the positive-energy vacuum with the 
smallest decay rate — still dominates the volume fraction. In particular, even if the dominant 
vacuum can yield classical transitions to another positive-energy vacuum, the volume fraction 
of the latter is still much less than unity. 

Generally speaking, the main effect of accounting for bubble collisions is to enhance the 
connectedness among vacua in the landscape. This is relevant to the so-called staggering issue, 
which concerns the competition between anthropic selection for small (magnitude) vacuum 
energies and cosmological selection for large volume fractions when attempting to explain 
the observed size of the cosmological constant, see for example [31, 32, [34[| . One might 
also take interest in the likelihood that our vacuum was created by a classical transition, as 
opposed to by semi-classical bubble formation. Our analysis indicates that this likelihood is 
determined by the detailed relationships among transition rates in the landscape. 

The remainder of this paper is organized as follows. In Section ^ we review the construc- 
tion of the rate equation, ignoring bubble collisions, while in Section |3| we include bubble 
collisions. We predominantly work in terms of a scale-factor-time foliation, though we briefly 
explain how to translate the results into those of a lightcone-time foliation. The phenome- 
nology of the rate equation is studied in Section ||. We begin with a review of a simple toy 
landscape, ignoring bubble collisions, and then we explore the toy landscape while including 
the effects of some representative classical transitions. We extend our results to a more gen- 
eral landscape in Section 4^, and discuss the likelihood of a classical transition in our past 
in Section 4.4. Finally, we draw our conclusions in Section |5l 



2. Rate equation without bubble collisions 

To begin, we reconstruct the standard rate equation, adopting the usual assumptions || [l9|] . 
In particular, we take the spacetime to be everywhere (3+l)-dimensional (the rate equation 
in a trans-dimensional multiverse is studied in |55|), we assume that all vacuum transitions 
occur via semi-classical bubble formation, we coarse-grain over the timescales of any transient 
cosmological evolution between epochs of vacuum-energy domination, we assume there are no 
vacua with precisely zero vacuum energy, 1 and we ignore bubble collisions. In the next section 

lr The scale- factor and lightcone-time cutoff measures both fail to regulate the volume in bubbles with precisely 
zero vacuum energy. (The author thanks Adam Brown for explaining the problem with the scale-factor cutoff 
measure.) Note that four-dimensional vacua with precisely zero vacuum energy might not exist. 
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we include bubble collisions, but the other assumptions are held throughout the analysis. 

The rate equation describes the volume fractions of vacua in terms of some global time 
variable t. To establish the global time foliation, we start with a large, spacelike hypersurface 
So on which we set t = 0. Note that it is not necessary for So to be a Cauchy surface: if So 
intersects an eternal worldline (that is, a worldline that never sees a vacuum energy density 
less than zero), then the rate equation describing the future evolution of So will possess an 
attractor solution that is independent of the detailed orientation and distribution of vacua on 
Sq. This implies that the volume fractions are independent of the choice of Sq. 

Suppose the total physical volume of So in vacuum i is V%. To construct the rate equation, 
we first compute the change in physical volume in vacuum i over a time interval At, i.e. 

m = Vi&At) - ^(E ) , (2.1) 

where E^ denotes the hypersurface of constant t = At. The coarse-grain approach explores 
the limit of large At to compute the various contributions to AVi, but then expands in 
At <C 1 to construct a differential equation. In the context of bubbles with positive vacuum 
energy densities — henceforth referred to as dS bubbles — this means that the rate equation 
ignores the transient cosmological evolution between epochs of vacuum energy domination. In 
the context of bubbles with negative vacuum energy densities — henceforth referred to as AdS 
bubbles — it means the rate equation ignores the cosmological evolution altogether. 2 Although 
the coarse-grain rate equation does not by itself provide an accurate assessment of the physical 
volume fractions in AdS vacua, it is convenient to track the volume fractions in these vacua 
anyway. We do this by simply conserving comoving volume during transitions to AdS vacua, 
and ignoring the subsequent evolution of the volume. 

2.1 Scale- factor time 



We first take t to be the scale-factor time 1 24 , 25 , 26, 27 



dt = Hdr , (2.2) 

where r is the proper time evaluated along a geodesic congruence orthogonal to So, and H is 
the local Hubble rate, in particular we can take H = (1/3) v?-^ in terms of the four- velocity 
field along the congruence. Although a precise definition of scale-factor time involves some 



subtleties in the treatment of locally contracting spacetime regions [ 36 , ^6], £7]] , these can be 
ignored in the coarse-grain analysis, which smears over such regions. 

Before proceeding, it is helpful to collect some facts about bubble formation. Consider a 
region in some dS vacuum i with cosmological constant Aj. On scales that are small compared 



2 This does not imply that the rate equation cannot be used to study cosmological evolution on timescales 
that are smaller than the coarse-graining, in dS or AdS vacua. The rate equation gives the bubble nucleation 
rate for all vacua that are reached primarily via semi-classical bubble nucleation starting from a dS vacuum. 
The bubble nucleation rate can then be combined with a more fine-grained analysis to study the dynamics 
in bubbles on smaller timescales. Some examples are reviewed in the appendix of |23|. 
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to the curvature radius, the line element can be written 

1 



ds' 



-dr z + 



4fl? 



2HiT ( dr 2 + r 2 ^ 



where 



(2.3) 



(2.4) 



cifi 2 is the line element on the unit 2-sphere, and the absolute value is inserted to establish a 
general definition for when we consider negative values of Aj. Now suppose an initially point- 
like bubble nucleates at time r = r nuc . To facilitate future reference we take the hypersurface 
T = r nuc to coincide with the aforementioned So, at least in the vicinity of the bubble. The 
bubble wall of a point-like bubble expands at the speed of light. Therefore, the comoving 
radius of the bubble at times r > r nuc is 



r w (r) = 2e 



-HiT n 



2e 



(2.5) 



The bubble expands so as to subtend a finite comoving volume (327r/3) e _3 ^ iTnuc in the limit 
r — > oo. Note that this comoving volume corresponds to a physical volume (4-7r/3) H^ 3 on the 
hypersurface So- Therefore, the loss of physical volume in vacuum i in the future evolution 
of So due to the nucleation of this bubble is equivalent (in the limit r — > co) to the loss of 
physical volume that results from simply ignoring the would-be future evolution of a physical 
volume (47r/3) H~ 3 on the hypersurface So- 

Now suppose that this bubble is a dS bubble. The coarse-grain analysis ignores dynamics 
on timescales smaller than the timescale of vacuum domination; therefore the line element in 
the bubble can generically be written 



ds 2 



dr 2 + —2 e 2H ^ T [di 2 + sinh 2 (£) dn 2 

3 



(2.6) 



where j labels the vacuum in the bubble, and it is implicit that we focus on bubble FRW 
times t S> HJ 1 . Importantly, surfaces of constant FRW time r in the bubble are not surfaces 
of constant scale-factor time t in the global foliation. In particular, it can be shown that 
the change in scale- factor time between the hypersurface So (corresponding to r = r nuc ) in 
vacuum i and a hypersurface of constant r in the bubble is |36| 



At 



Hjt + HHi/Hj) - ln(2) + £ + ^lnQ + ±e"« 



(2.7) 



Therefore, the induced metric on Sa< when it overlaps with the bubble is 



ds 2 



„2At 



H 2 



1 1 
- + -e 

2 2 



-4/3 



-2? 



+ H? \2 + 2 



-4/3 



e~ 2f sinh 2 (e)dfi 2 



(2- 
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The coarse-grain approximation explores the limit where the first term in brackets dominates 
over the second (because the coarse-grain approximation implies for example At > Hi/Hj). 
Accordingly, the physical three-volume of the intersection of £a* and the bubble is 

^sinh 2 (e) -jj Q + \e^ 2 e ~* = ^ . (2.9) 



Technically, the integrand in (2.9) is invalid at large values of £, for which the hypersurface 
^Ai explores times r < -^j" 1 an d the above results receive corrections. Nevertheless, for 
sufficiently large values of At the physical volume on S^t is dominated by regions for which 
r S> H~ , and the above approximations are accurate. 

Recall that we modeled the initial bubble as point-like. A realistic bubble has some non- 
zero initial radius, and the bubble wall has zero initial velocity. However, the bubble wall 
accelerates, its velocity approaching the speed of light. Therefore, in the limit of large At, 
our results for point-like initial bubbles coincide with the results for more realistic bubbles. 
An important exception to this rule is an "upward" transition from lower to higher vacuum 
energy, for which the bubble wall fills the horizon [||]. Lacking any clearer guidance, we simply 
use the above results for these vacua as well. 

We now return to the rate equation. Breaking the calculation of AVi up into parts, we 
first compute the change in physical volume in comoving regions that begin and remain in 
vacuum i over the interval At. This is given by the definition of scale-factor time: 

AVi = (e 3A * -l)Vi -> 3Vi At , (2.10) 

for any dS vacuum i. As mentioned above, we simply ignore the change in physical volume 
in comoving regions that begin and remain in the same AdS vacuum, as these vacua collapse 
into a big-crunch singularity on timescales that are small compared to the coarse-graining. 
Hence, for these vacua we take AVi = 0. 

Next, we compute the change in physical volume due to vacuum decay in comoving 
regions that begin in vacuum i. Referring to the above analysis, we note that when i is a dS 
vacuum then the effect of such decays is equivalent to removing physical volume (47r/3) H~ 3 
for each bubble nucleation, at the time of its nucleation. The number of such nucleations in 
an interval At is equal to the four-volume in vacuum i in that interval times the decay rate, 



AiVdccay = Yl T 3i V i AT ( A *) = Yl JT V i At ' 



Hi 

3 3 

where Tji denotes the transition rate from vacuum i to vacuum j, per unit physical three- 
volume per unit proper time, and we have used At = Hi At in vacuum i. The corresponding 
change in volume in vacuum i is therefore 

A ^ = -Ef§^ A *- ( 2 - 12 ) 
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This result assumes that i is a dS vacuum. Although AdS vacua can also decay |37|| , when 
they do so they transition exclusively to other AdS vacua, and therefore in the coarse-grain 
approach these transitions can be ignored (i.e. we set Tij = for AdS j). 

Finally, we compute the change in physical volume due to transitions to vacuum i from 
comoving regions that begin in some other vacuum. The number of such transitions is com- 
puted in analogy to ( [2,11 ), but now with reference to the transition rate from some (dS) 



vacuum j to vacuum i, summing over the relevant vacua j. In the limit of large At, each such 



transition generates a physical volume (2.9) in vacuum i, when i is a dS vacuum. This can 
be interpreted as the immediate creation of a physical volume (Att/3) H~ 3 , times a growth 
factor e 3A * which is already accounted for in ( |2.10| ). Thus we write 

A ^ = £||^ A *> ( 2 - 13 ) 

7 ■? 



for a dS vacuum i. For AdS vacua i, the analysis surrounding (2.9) no longer applies. Note 
however that the physical volume created in i due to transitions from a given dS vacuum 
j in ( 2.1 3| ) is precisely the same as the physical volume lost in j due to transitions to a 



given vacuum i m ( gTp (which requires exchanging the indices i o j). Since both of these 



expressions have the volume expansion factors stripped away, this equivalence expresses the 
conservation of comoving volume in vacuum transitions. Extending this principle to AdS 
vacua indicates that we should us ( |2.13| ) when describing the creation of volume in vacua i 
for dS and AdS vacua i. We make an error for AdS vacua since we ignore transitions from one 
AdS vacuum to another; however since these transition rates are exponentially suppressed 
next to the timescales of the AdS big crunches, this error is small. 
Combining results and taking the infinitesimal limit, we obtain 

^ = 3Vi + Y, K v y j ~ Yl K i iVi ' ( 2 - 14 ) 

j j 

where the first term is understood to apply only when i is a dS vacuum, and we have defined 
the dimensionless decay rates 

- 47rr u 

K ij = • (2-15) 

It is convenient to also define the volume fractions fi, where 



Since the dS vacua dominate the physical volume in the future evolution of Sq, we can write 



d Vj = %Y2j Vj, which means the rate equation can be written 



(if, 
dt 



^T = M ij f j , (2.17) 
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where the transition matrix is Mjj = Kij — 5{j K^i- The solution to ( 2.17 ) is |0J 

/<(*) oc /; 0) + Si e" 9 * + . . . (anti-de Sitter) 

(2.18) 

oc Si e" 9 * + . . . (de Sitter) , 

where the fjy are constants reflecting the initial conditions, q > is (minus) the smallest- 
magnitude eigenvalue of Mij, S{ is the corresponding eigenvector — called the dominant eigen- 
vector — and the ellipses denote terms that fall off faster than e~ qt . Although we have pre- 
sented the solutions for both dS and AdS vacua, as we have remarked the coarse-grain rate 
equation is does not reliably assess the volume fractions in the AdS vacua. 

2.2 Lightcone time 

The above analysis is sensitive to the choice of scale-factor time as the global time parameter 
t. Another popular choice is lightcone time ]15| , |i~6[ , defined according to 

t = -lln{V [l + (p))}, (2.19) 

where Vo[I + (p)] is the volume on So subtended by the subset of an initially uniform geodesic 
congruence orthogonal to So that intersects the causal future I + (p) of the point p at which 
the lightcone time is being evaluated. One can think of the geodesic congruence as a tool 
to project the asymptotic comoving volume of I + (p) back onto So- Thus, as one considers 
points p progressively further to the future of So, this projection covers a progressively smaller 
volume on So, and the lightcone time increases according to minus one third of the logarithm 
of this (another, equivalent definition is given in jig] ). 

Rather than repeat the full analysis of the previous subsection, we simply refer to the 



results of [18|. The crucial difference in comparison to scale-factor time is in the analogue of 
( [2.7D , i.e. the change in lightcone time between So and a hypersurface of fixed FRW time r 
in a bubble that nucleates on So- In the case of lightcone time, this is |l~8| 3 

- (gj + ifre-*) 4 - 
. (Hi + Hj)Hf 

Proceeding as in the previous section, the full consequence of this difference is to multiply 
fl2.9|) by a factor of Hf/Hj, which corresponds to multiplying the second term in ( p. 14 ) by a 



At = HjT - ln(2) + f + | lnQ + \e~^j 



(2.20) 



factor of Hf/Hf. Thus, we recover a rate equation of the form (2.14) if we work in terms of 



3 

the rescaled volumes V% = HfVi. Likewise we define 

7 _ Vi _ HfVj 

fi EjvrxsW (2 - 21) 

Since in the coarse-grain approach the various Hi are simply constants, this gives an equation 
of the form fl2.17p , from which the solutions can be read off. 

3 The analysis of [j^ does not keep track of all of the terms that we include here; for a more thorough analysis 
see That paper actually studies CAH time, but in the coarse-grain limit CAH time is proportional to 

the exponent of lightcone time, allowing us to convert from the results of that paper. 
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Figure 1: Left panel: toy spacetime diagram of a bubble collision producing a single domain wall. 
The bubble walls are represented by solid curves following null rays until they collide, producing a 
domain wall. Regions in vacua i, j, and k are labeled; the causal future of the collision is bounded by 
dotted lines and shaded gray. Right panel: the same as in the left panel, except the collision produces 
two domain walls (a classical transition) enclosing vacuum I. 



3. Including bubble collisions 

The calculation of Section ^ neglects the effects of bubble collisions, which we now address. 

Consider the collision between one bubble of vacuum j and another bubble of vacuum 
k, the two bubbles having nucleated in vacuum i (see in Figure ||). We assume that both 
j and k have smaller vacuum energies than i. When the bubbles collide, either the bubble 
walls annihilate (if j = k), or one or two domain walls form after the collision. In the case of 
one domain wall, the causal future of the collision contains both vacua j and k, their relative 
volume fraction determined by the trajectory of the domain wall. In the case of two domain 
walls, the domain walls contain some different vacuum t ^ j, k. If the vacuum energy of t is 
larger than the vacuum energies of j and k, the domain walls accelerate toward each other. 
Nevertheless, if the vacuum energy of I is less than the vacuum energy of i, the domain walls 
do not necessarily collide before future infinity. This is called a classical transition to vacuum 
£ 0, ||, H, H; in this case the causal future of the collision contains vacua j, k, and £, with 



the relative volume fractions determined by the trajectories of the domain walls. 

We use Xikji to denote volume fraction in vacuum I in the causal future of a collision 
between bubbles of vacua k and j which nucleated in vacuum i. As remarked above, for 
any given collision this volume fraction depends on the trajectories of domain walls. These 
depend on the model-dependent interaction potential governing the tunneling fields, as well 
as on the collision-dependent placement of the colliding bubbles. We are uninterested in these 
details, and instead simply take the Xikji as given, taking the given quantities to average over 
the relative placement of colliding bubbles (further details are presented below). 

3.1 Volume corrections from bubble collisions in scale- factor time 

Consider a dS bubble of vacuum j, which nucleates in some vacuum i at t = 0. According to 
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Q2.9D , the physical volume in this bubble on the constant scale-factor hypersurface Sa* is 

4tt 



V 



^3At 



(3.1) 



for large At. However, some fraction of this volume is not actually in vacuum j, because it 
resides in the causal future of collisions between the bubble of vacuum j and other bubbles. 
To account for this, we first compute the total volume in the causal futures of these collisions, 
so that we can subtract this volume of j from the rate equation. We then discuss the volume 
of each vacuum that should be put back into the rate equation so as to reflect the vacuum 
composition of the causal futures of these collisions. 

Part of the calculation is laid out nicely in |^] , and we begin by translating the relevant 
result into our notation. To do so, consider a 2-sphere of radius £ on a constant FRW time 
slice in the bubble geometry (|2.6| ). According to the symmetries of the collision, if the causal 
future of a collision intersects this 2-sphere, the intersection corresponds to a disk on the 2- 
sphere. This disk subtends a certain solid angle, and we are interested in the total solid angle 
subtended by summing over the causal futures of all bubble collisions. This is the calculation 
performed in |j{3] , and the resulting solid angle (in the limit of late FRW time slices) is 



4tt/(£) = 4vr 



1 



where 



n 



tot 



(0 = 4tt£ 



K ki 



+ In 



-JWO/4-7T 



S + 2^cosh(£) + l 



(3.2) 



(3.3) 



We 



where the sum runs over all vacua k into which the vacuum i can decay (with A& < A 
have used tan(T co ) = (Hi/Hj)tanh(Hjr/2) — > Hi/Hj in converting from the notation of 3£] 
in the context of the coarse-grain focus on the late-FRW-time limit in the bubble. 

To obtain the physical volume on E^t that intersects these causal futures, we integrate 
over 2-spheres of radii £, but using the solid angle (3^) instead of the usual Air, and using 
the induced metric (2.8) to switch from a constant-r to a constant-i hypersurface. This gives 

-2 



c^47rsinh 2 (£)#. 



-3 £ 3At 



1 1 

- + -e 

2 2 



-3 £ 3At 



F, 



(3.4) 



where the second expression defines F, which is a dimensionless function of R = Hi/Hj and 
Ki = Y2k K ki- The integral in F can be evaluated in terms of hyper geometric functions, but 
the result is not very illuminating. Instead we focus on the situation where R 3> 1, and we 
assume k» <C 1. Then we can approximate 

,-2 



F 



+ 



d^A-K sinh 2 (£) 

ln(fl) 

d(4n sinh 2 (£) 
d£ 4-7T e" ? 



1 1 

2 + 2 6 



1 



R 2 + 2i?cosh(£) + 1 



1 1 

- + -e 

2 2 



-t 



-3C 



1 - R~ Ki e 



(3.5) 



(3.6) 
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3 k 

where we have kept only the leading-order terms in Hi/Hj and Kij. Combining with ( |3.4| ), 
we find that the physical volume in the causal future of these bubble collisions is 

hT ( . :iA '>>,. (3.8) 



3H] Hi ^ 

Note that the volume subtracted from a bubble of vacuum j due to a collision with a 
bubble of vacuum k is not symmetric with respect to the indices j and k. In other words, 
when a bubble of vacuum j collides with a bubble of vacuum k, the causal future of the 
collision intersects a different "would-be" volume on SAt in the bubble of vacuum j, than the 
"would-be" volume on Xa* that it intersects in the bubble of vacuum k. There is no reason 
these two volumes should have been the same, because neither of them represents the actual 
geometry in the causal future of the collision; instead they represent the volume that would 
have been there had these bubbles not collided, and this counter- factual volume is different 
between the two bubbles because the two bubbles have different Hubble rates. 

After we subtract the volume ( |3.8| ) from a bubble of vacuum j, we must specify how much 
volume of each type of vacuum to put back in its place. We write the volume of vacuum t 
on the hypersurface £a* in the causal futures of the collisions between a bubble of vacuum j 
and other bubbles of vacuum k, for all k, as 

47T 



3Hg H t 



e 3At 



y] KkiXikji , (3.9) 



where i labels the vacuum in which the bubbles of vacua j and k nucleate. Of course, all of 
the microphysics is contained in the Xikjh which are also understood to reflect the average 
volume fractions after considering the various possible relative placements of the colliding 



bubbles that contribute to each term in (3.9). The other factors in ( |3,9| ) are chosen so that 
if collisions with bubbles of a given vacuum k always produce a domain wall that runs along 
the causal future of the collision into the bubble of vacuum k — that is, from the perspective 
of the bubble of vacuum j, collisions with bubbles of vacuum k have no effect on the interior 
of the bubble of vacuum j — then Xgkji = $ej (this returns all of the volume that had been 
removed from the bubble of vacuum j, but none of the volume that had been removed from 
the bubbles of vacuum k). 

The result ( p 7 ^ ) assumes that both % and j are dS vacua. Since we ignore transitions out 
of AdS vacua, this assumption about i is sufficient. On the other hand, it is important to keep 
track of the case where the vacuum j is AdS. This is because it is possible for the collision 
between two AdS bubbles to create a classical transition to a dS vacuum. Guided by the 



4 The vacuum reached by a classical transition must have a smaller vacuum energy than the vacuum in which 
the colliding bubbles nucleate. Therefore, two (AdS) bubbles that nucleate in an AdS vacuum cannot collide 
to produce a classical transition to a dS vacuum. This is one reason why it remains unimportant to track 
bubble nucleations in AdS vacua. 
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results of Section ||| in which conservation of comoving volume implied that after stripping 
away the volume expansion factors the same expressions could be used for transitions to dS 
and AdS vacua, we assume that fl3,8|) can be used to describe the volume in the causal future 
of bubble collisions for both dS and AdS vacua j, modulo the volume expansion factor e 3A *. 
Likewise, we use (|3.9| ) for both dS and AdS vacua £ (and j and A;). Although it is possible 
that using (|3.8|) for AdS vacua j introduces an error, technically this error can be removed 



by an appropriate choice for the factors \ekji m (3.9). Our analysis is unconcerned with such 
details; in our analysis it is only important whether \ikji is precisely zero or not. 

3.2 Modifying the rate equation in scale factor time 

We can now compute the change in volume in vacuum i, AVi, after a scale-factor time interval 
At, due to bubble collisions. We first consider the loss of volume in bubbles of vacuum type 
% due to collisions with other bubbles. This is 

AV * = ~ E jf. V i At ^Jf. K K = ~ E H '! h ^JT> V i At ■ ( 3 - 10 ) 

j,k ^ i J jfc i 

The first three terms in the first expression give the number of bubbles of vacuum i that 
nucleate in vacuum j in the time At, according to (2.11). The other terms in this expression 



give the volume in each such bubble that is in the causal future of collisions with bubbles of 
vacuum k, given by (|3.8| ), but where we have removed the expansion factor e 3A *, which in 
the coarse-grain analysis is accounted for by the factor of 3Vi in ( p. 14 ). 



Next we consider the change in volume in vacuum i, due the physical volume in vacuum 
i in the causal future of bubble collisions. This can be written 

AVi = E jt. V i At zh 2 H Khj Xlkej = E XiM i Kk i Ki i Jf2 V i At ■ ( 3 - n ) 

j,k,t i i i j,k,t i 

The first three terms in the first expression give the number of bubbles of type I that nucleate 
in vacuum j in the time At, according to fl2.ll ). The other terms in this expression give the 



volume in vacuum i in the causal futures of the collisions between each such bubble of vacuum 
i and bubbles of vacuum k, given by ( |3.9| ), but where we have removed the expansion factor 
e 3A *, which in the coarse-grain analysis is accounted for by the factor of 3Vi in ( 2.14 ). 
Putting everything together and taking the infinitesimal limit, we obtain 

-^=3Vi + Yl K ii V J ~ E K H V i + E X >l''J ,; lo f ''jJpi V j - E ' , ' ; .< Ki '.'777 r ; • ( 3 - 12 ) 

3 j j,k,l 1 j,k 1 

In terms of the volume fractions fi defined in ( |2.16| ), we have 

f = M tJ f 3 , (3.13) 

where now 

Mij = Kij — 5{j Kki + liktj K kj K lj i 

(3.14) 
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where for later convenience we have defined 

_ Hf 

liklj = Jj2 V^iklj ~ ■ (3.15) 

Note that conservation of comoving volume implies that the sum over rows in any column of 
Mij must be zero. This implies a useful constraint on the ji^ij, namely 

^2lik£j = 0- (3.16) 



3.3 Lightcone time 

The previous two subsections worked in terms of scale-factor time. However, the only place 
where the actual definition of the global time parameter enters is in the induced metric used 
in ( |3.4p (and results, such as ([3.1|), taken from Section |2|). Meanwhile, in Section |2.2| we 
found that the only effect of using lightcone time instead of scale-factor time corresponded to 



multiplying the second term in (2.14|) by H? /Hf. The same applies here with respect to the 



second term in ( |3.12 ), and because of the modification to ( |3.4D we multiply the fourth and 
fifth terms in ( |3.12j ) by Hj /Hf as well. It is easily checked that if we define /, and Vi as in 



( |2.21 ) and above it, we obtain 

dt 



Mijfj , (3.17) 



with Mij given by ( 3T3 ). Therefore, given a solution in terms of scale- factor time, we can 



again read off the solution in terms of lightcone time. 



4. Phenomenology 



To develop intuition for the phenomenology of the rate equation, we study a simple toy model 
of the landscape [19, [E], 26]. The model can be represented by the diagram 



1 «-> 2 3 

I ; ; . (4.i) 

4 5 6 

The numbers label different vacua, while the arrows indicate the direct transitions that are 
allowed among the vacua. (We use the phrase "direct transition" to designate semi-classical 
bubble formation via quantum tunneling through one potential barrier, which we take to be 
the dominant form of vacuum transition, aside from perhaps classical transitions, which are 
discussed below.) For concreteness, we assume the vacuum energies of this model obey 

A 4 , A 5 , A 6 < < Ai, A 3 < A 2 • (4.2) 

That is, the labels "1," "2," and "3" designate dS vacua, while the other labels designate AdS 
vacua. We also assume that "1" corresponds to the dS vacuum with the smallest decay rate. 
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4.1 Toy model without bubble collisions 

We first ignore bubble collisions. Then this model is similar to models studied elsewhere in 



the literature [19, 40, 26]. Seeking the late-time attractor solution — that is, the dominant 



eigenvector — we insert the ansatz /j,/, = s^e qt into the rate equation, obtaining 

— qs\ = K12S2 - K\si (4.3) 
-qs 2 = K21S1 + K23S3 - K2S2 (4.4) 
-QS3 = K32S2 ~ k 3 s 3 , (4.5) 

where we focus on the volume fractions of the dS vacua, and «j = ^ ■ Kji. This system of 
equations can be solved algebraically, but the solution is complicated and not very enlight- 
ening. Since K21 an d K23 are transition rates from lower to higher vacuum energy, they are 
expected to be exponentially suppressed relative to the other transition rates. We therefore 
solve the above system of equations by expanding in K21 and K23, treating them as similar in 
order for the purpose of the expansion. 

At zeroth order the only non-zero dS component of Sj is s±, which we can set to unity. 
The corresponding value of q is 

q = K41 + K51 = «i , (4.6) 

where the second equality implicitly evaluates K\ at zeroth order, a notational shorthand that 
we also employ below to help simplify expressions. The other dS components of s, become 
relevant at first order, and are given by 

«21 K21K32 , , 

S2 = , s 3 = . 4.7) 

The zeroth-order solutions for s± and q receive first-order corrections, but these are subdomi- 
nant and so we ignore them. We see that the dS components of s» are dominated by "1" — the 
dS vacuum with the smallest decay rate — which is accordingly called the dominant vacuum 
(this assumes that the decay rates are not tuned so as to make q very near to K2 or k%, in 
which case there could be a set of degenerate dominant vacua ]2f| ) . Intuitively, the dominant 
vacuum dominates the asymptotic volume of spacetime due to its small decay rate, and the 
volume fractions of all of the other vacua reflect in part the relative likelihood of transitioning 
to them from the dominant vacuum. In particular, S2 carries a factor of K21, reflecting the 
transition from "1" to "2," while S3 carries the factor /t2i K 32j reflecting transitions first from 
"1" to "2," then from "2" to "3." These factors also contain total decay rates in the denomi- 
nator, so they only correspond to suppression factors when the corresponding transition rate 
is relatively suppressed. This applies to K21 in the first-order terms in (fO|), and it applies 
to contributions to the s, that arise due to transitions through "3," which are suppressed 
relative to the above contributions. 

4.2 Effects of bubble collisions 

The above discussion summarizes previous work. Our goal is to understand the implications 
of allowing for bubble collisions. Collisions that result in a single domain wall, and therefore 
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merely shift the volume fractions of the vacua involved in the collision by order- unity factors, 
are not consequential in the context of the qualitative dynamics described above. However, 
the possibility of classical transitions might change the above conclusions in more dramatic 
ways. To explore this, we first assume that bubbles of vacua "4" can collide and produce a 
classical transition to vacuum "3." We denote this with the diagram 

4)(4-»3. (4.8) 



Note that this transition requires that A3 < Ai 11, 2£]. Although the landscape ([4.1[) per- 
mits other bubble collisions, for the moment we assume that these do not produce classical 
transitions so that we can ignore them. Then the full set of non-zero %jkt can be taken to be 

ff2 tt2 

1 1 , . 

74441 = -p> 73441 = ' 

The first term accounts for the subtraction of volume in the causal future of the collisions 
that would be in vacuum "4" but for the collision, while the second term accounts for the 
restoration of this volume but in vacuum "3," a consequence of the classical transition. In a 
realistic model we expect some of the causal future of the collision to be in "4," which means 
that 74441 should be less negative and 73441 should be smaller. However, these effects change 
our results by factors that are negligible next to the exponential staggering of decay rates, 
and so for simplicity we ignore them. 

Including these collisions and their effects, the rate equation gives 

-qsi = ki 2 s 2 - K4S1 (4-10) 
-qs 2 = K21S1 + K23S3 - k 2 s 2 (4.11) 
-QS3 = k 32 s 2 - K3S3 + 73441 K41S1 , (4.12) 

where again we focus on the dS components of Sj. As before, the solution is more transparent 
if we expand it in terms of K21 and K23. At zeroth order, we find 

S1 = ( 1 + m^kY 1 , S3 = i^iA , ( 4.i3) 

V «3 - q J K3-q 

where we have normalized the zeroth-order solution so that Y2i s i = 1> where the sum runs 
over dS vacua i, and q is unchanged from before: q = m. The leading-order contribution to 
s 2 still appears at first order in the expansion, 

S 2 = ^i£L + ^£3_ (414) 
K 2 -q K 2 -q 

where s± and S3 refer to the zeroth-order quantities. Evidently, the hierarchical structure 
of the dominant eigenvector is qualitatively changed relative to before. In particular, while 
the dominant vacuum "1" still dominates the dS components of Sj, the component S3 now 
contains a product of two downward transition rates from "1," as opposed to an upward 
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transition rate. Note that since "1" is by definition the dS state with the smallest decay rate, 
we still expect S3 to be exponentially suppressed relative to s±. Nevertheless, the size of the 
suppression is dramatically reduced. Meanwhile, the volume fraction of "2" is still suppressed 
by an upward transition rate, but this can come from "1" or from "3." 

Although the results are changed from when we ignored bubble collisions, our intuition 
for the solution is maintained. In particular, the dS vacuum with the smallest decay rate 
dominates the asymptotic volume fraction, and the volume fractions of the other dS vacua 
reflect the relative likelihood of transitioning to them from the dominant vacuum. In the case 
of "3," this likelihood is enhanced by the possibility for classical transitions. This in turn 
enhances the volume fraction of "2," insofar as it can be reached by transitions from "3." 



The classical transition represented by (4.8) is not the only possibility in the landscape 



(4.1). Another possibility occurs when bubbles collide in vacuum "2." We now focus on the 
possibility that collisions between bubbles of vacua "5" can result in a classical transition to 
"3," represented by the diagram 

5)(5 — »> 3, (4.15) 
and ignore all other possible collisions. Then the full set of non-zero "fijkt can be taken to be 

H 2 H 2 

75552 = -77I ' 73552 = 77§ , (4.16) 

^3 -"3 



where the discussion below (4.G) applies here as well. The resulting rate equation gives 



-qsx = Ki 2 s 2 - K4S1 (4-17) 
-qs 2 = «2isi + K23S3 - K2S2 (4.18) 
-QS3 = k 32 s 2 - K3S3 + 73552K52S2 • (4.19) 

We again expand in terms of n 2 i and n 2 3. The zeroth order dS components are the same as 
without classical transitions, namely the only non-zero component is s\, which we can set to 
unity. The other dS components of become relevant at first order, and are given by 

"21 («32 + 73552«52) K 2i lA on x 

S2 = , s 3 = — — . (4.20 

K 2 -q {K 2 -q){K3-q) 

Note that s 2 and S3 are still suppressed by an upward transition rate out of "1." Therefore, 
the effects of classical transitions do not change the qualitative expectations for the dominant 
eigenvector described above. On the other hand, it is possible for the component S3 to be 
enhanced relative to before, depending on the relative size of 73552*52 and K32. If K32 were 
zero then the connectedness of the landscape would be greatly enhanced. 

4.3 Generalization of the toy model 

The toy landscape model Q4.1|) is simple, but the intuition developed above extends to more 



general landscapes. Indeed, it is possible to compute the components of the dominant eigen- 
vector Si for any number of vacua and for any set of transitions among them, by expanding 
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in the off-diagonal elements of the transition matrix My p3| , [iT|| . As before, we assume the 
dS vacuum with the smallest decay rate is unique and denote it as "1." The (unnormalized) 
dS components of s, are then, to leading order, 

si = y _^ !£L _ x _ _ _ x K y l^ lx _ x S i i ] (421) 

{Pa} K{ ~ q ^"9 q 

where q is the decay rate of "1" and the sum covers the sequences of transitions that connect 
"1" to i using the fewest number of upward transitions. In the second expression we have 
exploited the fact that q is the smallest decay rate to write K{ — q ~ Kj, which allows for the 
terms to be rearranged so as to express factors as branching ratios. 

We expect the dS vacuum with the smallest decay rate to also possess a small vacuum 
energy (for some discussion of this matter see [pq| ). This ensures that transitions from "1" 
to other dS states are either transitions to larger vacuum energy or transitions that involve a 
very small change in vacuum energy compared to the expected size of the barrier. The rates 
associated with such transitions are exponentially suppressed relative to the rates associated 
with transitions to AdS vacua with much larger magnitude vacuum energies M. Therefore, 



we expect the right-most term in (4.21) to be exponentially small. This maintains the notion 
of "1" as the state that dominates the volume fraction, i.e. the dominant vacuum. 

It is straightforward to generalize the above result to include the effects of bubble col- 
lisions. In particular, the effect of bubble collisions is simply to add the (small) quantity 
Ylk ilikij^kj^ij to the quantity «y corresponding to each off-diagonal perturbative element 
of the transition matrix My. We can then read off the results from the analysis of [j33| , 41]: 



Si 



E K «Pl+Sj,fc7«jfcpi K Jpi K fcpi Kp n l + Ylj, k 7p„jfcl K jl K fcl / A nn\ 
x ... x (4.22) 



{M 



Ki q Kp n q 



Eq K ipi + Ylj,k Jijkp 1 Kjp 1 Kkp 1 K p n l + J2j,k lp n jkl K jl K kl 

— - X ... X - . 

r i K « K pi y 

{Pa} 

Of course, it is not necessary for both terms in a given factor to be non-zero, and in fact 
the consequences of bubble collisions are most significant when the first term in one of these 
factors is zero, as with the classical transitions described above. 

We find that classical transitions have enhanced the connectedness of the landscape, by 
which we mean in an expansion in terms of the number of upward transition rates -/V up , the 
number of non-zero components of Si at a given order in A^ up is increased due to classical 
transitions. This is relevant to the staggering issue, which concerns the competition between 
anthropic selection for small (magnitude) vacuum energies and cosmological selection for large 
volume fractions when attempting to explain the observed size of the cosmological constant 
[ |3l| , 32, 33, 34]. In particular, one aspect of the landscape that is important to the staggering 



issue is the rate at which the number of non-zero components of Sj increases as one considers 
larger values of A^ up the expansion in upward tunneling rates. Classical transitions allow 
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this rate to be significantly larger than otherwise. Many additional considerations enter the 
analysis, however, so further exploration of this issue is beyond the scope of this paper. 

If an anthropic vacuum could be reached via classical transitions caused by AdS bub- 
ble collisions in the dominant vacuum, it would imply dramatic consequences for landscape 
phenomenology. The prevalence of this vacuum would be exponentially larger than the preva- 
lences of anthropic vacua that arise only after upward or lateral transitions from "1," thereby 
allowing for sharp cosmological predictions (this is essentially an extreme case of the stagger- 
ing issue mentioned above). However, we consider this possibility very unlikely. A large land- 
scape of A diverse vacua seems most natural in the context of a roughly log(A)-dimensional 
configuration space |42"| , 43, 44]. Thus, we assume that each vacuum can directly transition to 
roughly log(A) nearby states. Classical transitions are permitted based on special relation- 
ships among local minima in the vacuum configuration space, such that each unique pair of 
vacua involved in a collision can classically transition to at most one unique state piEasa. 
Hence, for each vacuum there are at most roughly log 2 (A) distinct states that can be reached 
by a classical transition. In the context of string theory, one usually takes log(A) ~ 0(1000) 
or so [42, 43]. While A is then more than large enough to (for instance) fine-tune the cosmo- 
logical constant to anthropically acceptable levels |4q| , log 2 (A) is far from sufficient. Indeed, 
it is not clear that the roughly log 2 (A) states (optimistically) accessible to "1" via classical 
transitions is sufficient to contain a dS state at all, given the constraint that "1" can only 
classically transition to states i that satisfy Aj < Ai. Roughly speaking, we expect this pos- 
sibility to be unlikely if the vacuum energy of "1" is such that GAi < l/log 2 (A). While the 
vacuum energy of the dominant vacuum "1" is unknown, there are indications that it easily 
satisfies this inequality for log(A) ~ 0(1000) (see again pGfl). 



4.4 Probability to observe a classical transition 

So far our discussion has focused on global volume fractions. On the other hand, if the initial 
conditions for the phase of slow-roll inflation believed to describe our past were established by 
a classical transition — as opposed to by the formation of a single bubble — then this could have 
consequences for observational cosmology, at least if slow-roll inflation did not last too long. 



Indeed, the collision geometry retains an SO(2,l) symmetry [46], and so the phenomenological 



consequences should resemble those of the "anisotropic bubble nucleation" scenario studied 



in [47, 48]. We leave further exploration of these signatures to future work, and here focus 
on assessing the relative likelihood of this possibility. 

One way to approach this problem is to compare the rate of classical transitions to our 
vacuum to the rate of bubble nucleations of our vacuum. The results of Section || allow us to 
go a step further, since there we have computed the physical volume on a fixed scale-factor 
time hypersurface in the causal futures of these events. In short, these are given respectively 
by the fourth and second terms on the right-hand-side of the physical-volume rate equation, 



( 3.12| ). Using Vi oc Sj e^ 3 we write the ratio of the volume in i in the causal future of a 
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classical transition to the volume in i in the causal future of a bubble nucleation 

R = = . (4.24) 

l^j K ij S j 

Evidently, R depends on the detailed microphysics of the landscape. In particular, while the 
sums over vacua j in the numerator and denominator both cover all dS vacua j, for almost 
all of these vacua the transition rates mj are zero, while the vacua j for which the rates are 
non-zero are in general different in the numerator and denominator. Since the components 
of the dominant eigenvector Sj are exponentially staggered, R could be enormous or zero 
depending on the detailed structure of a large set of transition rates in the landscape (zero 
corresponds to the case where classical transitions cannot create vacuum i). 

It is amusing to consider a highly idealized situation. Suppose that classical transitions 
to our vacuum can only be caused by bubble collisions in one type of vacuum, and suppose 
that bubbles of our vacuum can only nucleate in one type of vacuum, and suppose that the 
Hubble rates and the relevant decay rates and the components of the dominant eigenvector 
for these vacua are the same: Hq, T, and so respectively. Finally, suppose that the classical 
transition converts the entire causal future of the collision into vacuum i. In this situation, 
( |3.16| ) and ( [3.15 ) imply that Ylki^iktj = 1- Putting all of this together, we find 

This expression is familiar from the bubble-collision literature, and gives the typical number 
of bubble collisions in an observer's past lightcone (see for example |^]). In this context, 
Hi is the Hubble rate during slow-roll inflation in the observer's bubble, while Hq and T are 
respectively the Hubble rate and vacuum decay rate outside of the bubble. The congruence 
between these two results is not surprising: the idealized situation that we have established 
above is equivalent to comparing the volume in the causal future of bubble collisions to the 
volume not, assuming that the volume in the causal future of two collisions is double-counted, 
etc. This is equal to the number of bubble collisions in an observer's past lightcone. The 
appearance of the inflationary Hubble rate Hi in the bubble-collision result, compared to 
what one would might suggest should be the Hubble rate associated with late-time cosmo- 
logical constant domination in the relative probability ( 4.25| ), is due to the former focusing 
on observers like us who arise at finite time while the later (interpreted as above) implicitly 
focusing on hypothetical observers at future infinity. 



5. Conclusions 

We have updated the rate equation describing the volume fractions of vacua in the multiverse 
to account for bubble collisions. The rate equation refers to a global time foliation, and we 
have presented our analysis in terms of a scale-factor-time foliation, also providing a couple 
of brief remarks on how to translate the results into those of a lightcone-time foliation. 
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As in the case where bubble collisions are ignored, the intuition for the asymptotic attrac- 
tor solution to the rate equation revolves around the dS vacuum with the smallest decay rate. 
Owing to its stability, this vacuum dominates the volume fraction of the multiverse — hence 
it is called the dominant vacuum — and the relative volume fractions of any other vacua are 
understood in terms of the relative likelihoods of the sequences of transitions that connect 
the dominant vacuum to them. The most significant effect of bubble collisions stems from 
the possibility for classical transitions to vacua that are otherwise not involved in the colli- 
sion. These transitions should be included along with bubble nucleations in the sequences of 
transitions referred to above. Thus, the main effect of accounting for bubble collisions is to 
enhance the connectedness among vacua in the landscape. 

The global volume fractions of dS vacua play an important role in determining whether 



a given measure has a Boltzmann brain problem [49, 5C, 51, [5^, 53]. Insofar as the effects of 
bubble collisions do not change the qualitative hierarchical structure among the components 
of the dominant eigenvalue, previous conclusions on this issue are maintained (namely, that 
Boltzmann brain domination can be avoided with plausible assumptions about the landscape) 
[ f2g| , |36[ ]. On the other hand, if the dominant vacuum yields classical transitions that create 
a dS vacuum, then this modifies the qualitative hierarchical structure of the dominant eigen- 
value, since the volume fractions of such vacua are not suppressed by an upward transition 
rate. We have argued that such transitions are unlikely to be possible, due to the expected 
smallness of the vacuum energy of the dominant vacuum. Nevertheless, even if dS vacua could 



be created in this way, they could only change the conclusions of [26, 36] if it were capable 
of supporting Boltzmann brains. We have argued that the number of vacua accessible to the 
dominant vacuum via classical transitions is at most roughly log 2 (N), where N is the number 
of vacua in the landscape. This implies that dS vacua created in this way are unlikely to have 
vacuum energies well below 1/G 2 log 2 (N). Meanwhile, [26] found that the largest vacuum 



energies capable of supporting Boltzmann brains must be well below (1/G 2 ) 10 -16 . Setting 
aside the other special circumstances necessary to create Boltzmann brains, we conclude that 
bubble collisions do not introduce a Boltzmann brain problem granted plausible assumptions 
about the landscape. 
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